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Abstract 

A simple unifying mass matrix is presented for the three active and one sterile 
neutrinos Ve, v^, v-r, and i^s, using an extension of the radiative mechanism proposed 
some time ago by Zee. The total neutrino-oscillation data are explained by the scheme 
T^e ^ T^s (solar), <-> Vr (atmospheric) and fg ^ (LSND). We obtain the inter- 
esting approximate relationship (Am^)atm — 2[(Am^)soiar(^"i^)LSND]^''^ which is well 
satisfied by the data. 



Three neutrinos, each associated with a charged lepton (e, /i, r), are now known. The 
invisible width of the Z boson, coming from the decay Z z/z/, is also consistent fll] with 
exactly three such neutrinos. This means that if there is a fourth neutrino, either it has to 
be very heavy (with mass greater than or it does not couple to Z. In particular, 

if it is light, then it must not have any electroweak gauge interactions. Such an object is 
often referred to as a "sterile" neutrino. The reason that this may be a necessary part of 
our understanding of particle physics is that there are at present three classes of neutrino 
experiments]^, 0, |^ which show evidence of neutrino oscillations with three very different 
Am^'s, i.e. differences of mass-squares. If all three interpretations are correct, then we need 
four light neutrinos. (A possible but rather extreme three- neutrino scenario is to have 
large anomalous z/T--quark interactions.) It is thus of theoretical interest to find a natural 
mechanism which explains the masses and mixings of these four neutrinos in the present 
experimental context. 

A specific model for a 4 x 4 neutrino mass matrix was proposed already some time 
ago.. The form of this matrix agrees with subsequent purely phenomenological analyses]^, § 
of all neutrino-oscillation data. Our present study concerns the possibility that all neutrino 
masses are zero at tree level, but are generated radiatively at one-loop to match the pattern 
in [6], using a mechanism first proposed by ZeelQ. We extend previous work|TO|, pT| on this 
topic to include a sterile neutrino with the help of an extra U(l) gauge symmetry [P^]. 
The resulting mass eigenvalues lead to the approximate relationship 

LSND \^) 

which is well satisfied by the data. 

Our model extends the standard electroweak gauge model to include three singlet fermion 
fields VsLi Nfi, and Sr, as well as 3 singlet scalar fields xt: xt i There are also two 

scalar doublets (05^,0°) and (02^,02); where only one is needed in the minimal standard 
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Table 1: List of fermion and scalar fields in our model. 



model. To obtain radiative masses for the three doublet neutrinos, just (</>]^2, 2) Xi 



are enough [TTJ. The more difficult task is to include the singlet neutrino into a 4 x 4 
radiative mass matrix of the same form. A natural way that this may come about is to 
have an extra gauge symmetry U{1)' for the fields Vsl-, xt ■> ^^^d X2 which is broken at a 
higher (~ TeV) scale. The axial- vector anomaly, generated by Vsl-, is cancelled by Nr which 
transforms as z/^l under U{1)'. We also add Sr which is trivial under U{1)'. A large mass 
for Nr is then ensured through the Yukawa interaction SrNr^x^ since {x%) >~ 1 TeV. The 
particle content of the model is summarized in Table 1. 

We have an unbroken discrete Z2 symmetry, namely L-parity, to distinguish between two 
classes of fermions. The leptons now have odd L-parity, replacing the usual additive lepton 
number. This allows the four neutrinos to acquire Majorana masses. However, tree-level 
neutrino masses are forbidden by the assumed particle content of our model, even after the 



spontaneous breaking of the gauge symmetry. Note that does not get a Majorana mass 
because of U{1)'] it also does not get a Dirac mass by pairing up with Nr or Sr because 
of L-parity. More specifically, consider the following interaction Lagrangian density of the 
fields shown in Table 1. 



i,j i i 

+l^{<Pt<p2 - 0?02 )xr + /^'X^X2 X2 + h'NRSnxT + h.c, (2) 



where we have used the notation ipiQ = ipi'Cj fo^^ two fermion fields and Evidently, 
fij is antisymmetric in its generation indices. We have assumed in (2) that (^2)^2) do 
not couple to leptons. This is easily achieved by a separate discrete Z2 symmetry which is 
explicitly broken, but only by soft terms such as + <t>\*<P\ + h.c. in the Higgs potential, 

as in the minimal supersymmetric standard model, for example. As shown below, the above 
interactions induce a radiative neutrino mass matrix for i/g, i^^u, Vr-i ^-nd Vg of the form 



Q a h d 

a c e 

b c f 

d e f 



(3) 



which generalizes the 3x3 matrix of the Zee model [9] by including a fourth row and column. 

In Fig. 1 we show the one-loop diagram linking Ui and Vj which contributes to the cor- 
responding entry in M.^- This is of course identical to that of Ref. [9] and [11]. Note that 
i j necessarily, hence only off-diagonal entries can be nonzero. Since hi — mij{4>1), 
obtain 



t I 2 2\ ( l^'^2\ 7-1/2 2 \ 



(4) 
(5) 
(6) 



where vi^2 = (0i2)) the function F is given by 



F{ml,ml) 



1 1 m? 
IGtt^ mf — TO2 ml 



(7) 



In Fig. 2 we show the analogous one-loop diagram linking i/^ to i/g- We find 



fl u 



(8) 
(9) 

(10) 



where u = (x2)- t^^® following, we will assume that fgUie is negligible. Moreover, while u 
is expected to be large compared to vi^2, that can be compensated by m^^ being larger than 
m-^^ or 777.0^. Thus d, e, / need not be larger in magnitude than a, b, c. In any case, we have 
the important relationship 

(„) 

c V rn^J fru 
where ml in Eq. (5) has been neglected. 



We make the same observation as in Refs. [9] and [11] that h and c are hkely to be the 
dominant entries of A^jy because they are proportional to m^. This means that Vt combines 
with a linear combination of Ve and to form a pseudo-Dirac particle. Let us also assume 
that |/er| << |//ir|) SO that |6| << |c|. Then the 2x2 submatrix spanning and Ug is given 

by 

— 2a6/c d — be/c — af/c 
d — be/c — af /c —2ef/c 

-2ab/c -ffeJUr - {be/c){ml/ml) 

-ffe,/Ur - (be/c) {mil ml) -2ef/c 

where we have used Eq. (11) and the fact that \a/c\ « {fen/ furl- Hence 



(12) 



and for rrii,^ « niu^, the Ue — Vs mixing is {cfen/efur + 6m^//m^)/2. This is assumed to be 
small, so as to satisfy the solar neutrino data. We now have 

(Am2),oia. ^ 4^. (14) 

Since M.^ has zero trace, it can easily be shown that the leading expressions for its 
eigenvalues are given by 

ah ah ef ah ef ^ef ,^ ^, 

-2—, c+ — + -c+ — + -2^. 15 

c c c c c c 

Hence the mass-squared difference between the two Majorana components of the pseudo- 
Dirac neutrino with mass c is 

Am^ = A{ah + ef) ~ 4e/ ~ (Am^jatm- (16) 

Since this is for a i/^ — i/,- mixing of 45°, we have taken it to explain the atmospheric neutrino 
data. Finally, the LSND data involve the mixing of and i^^, hence 

(Am2)LSND = c^ (17) 

with mixing given by h/c. Combining Eqs. (14), (16), and (17), we obtain Eq. (1), as claimed. 

Current neutrino-oscillation data are consistent with (Am^)LSND ~ 1 eV^ and (Am^)soiar 
~ 6 X 10"^ eV^. In that case, (Am^)atm is predicted by Eq. (1) to be about 5 x 10"^ eV^, 
which is supported by the most recent data from Super-Kamiokande. In our model, i/^ and 
Vt have the same mass c ~ 1 eV and they mix maximally. Let 6 ~ 0.04 eV, then the i^^ — Ve 
mixing parameter (sin^ 26')lsnd is 'ih'^/c? ~ 6 x 10~^, in good agreement with data. For 
Ve — oscillations, we let 



SO that (sin^ 2^)soiar is also about 6 x 10~^, again in good agreement with data. More 
specifically, we can let e ~ 0.12 eV and / ~ 0.01 eV, then rriv^ ~ '^&f /c ~ 2.4 x 10~^ eV. 



Furthermore from Eq. (18), fen/ f^r is now about 0.008 and from Eqs. (4) and (6), a ~ 3xl0~^ 
eV, hence m^^ ~ 2ab/c ~ 2 x 10~^ eV, justifying our assumption that m^^ « my^. We have 
thus a completely successful phenomenological picture of neutrino oscillations. 

The model of Ref.[ll] differs from ours in that Vg is assumed there to acquire a tree- level 
mass which is just slightly bigger than the radiative mass of v^- [This is of course rather ad 
hoc, but it is necessary to satisfy solar data.] Let us compare its consequences with those 
of ours. In the former, the parameter a is forced to be large in magnitude because 4a6 is 
identified there with (Am^)atm, resulting in |/e^| ~ iZ/irl- This condition is subject to severe 
phenomenological constraints because /e^ contributes to // decay. In fact, in that scenario, 
flii ^ fir ^ 10~"^Gf- (cos^ 0Mf^ + sin^^Ma"^)"^ where Mi,2 are the physical charged 
Higgs masses and is their mixing angle. In our model, because of Eq. (16), a can be and 
is very small, hence \feij\ « I/mtI, so that our is not constrained to be small. 

We note also that the form of Eq. (3) for the neutrino mass matrix with c as the dominant 
entry is not sufficient by itself to have the correct Ve—fs submatrix needed to explain the solar 
data. Without Eq. (11), which is an automatic consequence of our model, that submatrix 
would have dominant off-diagonal terms, i. e. 



which would make and pseudo-Dirac partners with the requisite mixing of 45° in conflict 
with solar neutrino data. 

A third point concerns the fermion singlets Nr and Sr. They have even L-parity, which 
is unbroken in our model, hence they do not mix into the lepton sector. Both of them are 
massive, because the terms SrSr and NrSr')^ -\-h.c. are allowed in the Lagrangian density. 
The scale of t/(l)'-breaking, i.e. (X2) can be taken beyond 1 TeV, thereby pushing up these 
masses. It is to be noted that the off-diagonal terms in the Z — Z' mass matrix are prohibited 
due to the absence of appropriate Higgs fields in the present model. We also assume that 




(19) 
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the kinetic mixing between the U{1)y and gauge bosons is neghgible. Hence our Z' 

couples at the tree level only to Usl, Nr, xt ^"^^ X2- Thus present experimental bounds [14] 
on a possible Z' with standard-model-like couplings do not apply. However, because z/^ mixes 
with Vf. radiatively, Z' develops a small coupling to z/g. To avoid any possible conflict with 
nucleosynthesis or current electroweak phenomenology, we assume Mz ~ 1 TeV or greater, 
which is of course natural since we already assumed {X2) ~ 1 TeV or greater. 

The charged scalar Xi contributes to the standard-model effective charged- current inter- 
action due to the presence of the fijivuliL — hL^jL)Xi term in the Lagrangian density. The 
corresponding effects on processes, such as electron-neutrino scattering, are experimentally 
severely constrained. They give rise to the constraint /g^/M^ < O.OSGGp [11], where M is 
the mass of the charged scalar mediating the process. Since we have \fe^\ « |/^(t|, this is 
no problem for us. The proposed hierarchical relation |/e^| << is also consistent with 
the constraint from the branching ratio of the decay fi'UnUr being (17.35 ± 0.10)% 



| 1^ since the latter only requires f^^/M < 0.13 Gp- 

In summary, we have demonstrated that the present results of solar, atmospheric as well 
as LSND experiments can be explained with three electroweak-active neutrinos and a sterile 
one with a minimal extension of the standard SU{2)i x f/(l)y electroweak gauge model. 
The extra U{1)' gauge and Z2 discrete symmetries are needed to avoid tree-level Majorana 
or Dirac mass terms. All neutrino masses are radiatively generated in one loop by an 
extension of the Zee model. Our proposal results in an interesting relationship (Am^)atm — 

1 /2 

2[(Am^)soiar(Am^)LSND] which is well satisfied by the present experimental data and will 
be critically tested with more accurate data forthcoming in the near future. 
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FIG. 1. One loop radiative z/j — = e,fi,T) mass due to charged Higgs exchange. 
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FIG. 2. One loop radiative z/j — Ug {i = e, /x, r) mass due to charged Higgs exchange. 
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